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Abstract 

The explosive crystallization of germanium ultradisperse amorphous films is studied 
experimentally. We show that crystallization may be initiated by local heating at the 
small film thickness but it realizes spontaneously at the large ones. The fractal pattern 
of the crystallized phase is discovered that is inherent in the phenomena of diffusion 
limited aggregation. It is shown that in contrast to the ordinary crystallization mode 
the explosive one is connected with the instability which is caused by the self-heating. 
A transition from the first mechanism to the second one is modelled by Lorenz system. 
The process of explosive crystallization is represented on the basis of the self-organized 
criticality conception. The front movement is described as the effective diffusion in the 
ultrametric space of hierarchically subordinated avalanches, corresponding to the explosive 
crystallization of elementary volumes of ultradisperse powder. The expressions for the 
stationary crystallization heat distribution and the steady-state heat current are obtained. 
The heat needed for initiation of the explosive crystallization is obtained as a function of 
the thermometric conductivity. The time dependence of the spontaneous crystallization 



probability in a thin films is examined. 
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1 Introduction 

The metastable amorphous films are obtained usually by quenching of a melt or by steam 
condensation on a cold substrate |]J -[[§]. Experiments show a vast variety and complicated 
character of the subsequent transitions into the stable crystalline state If the 

film thickness is so small that the crystallization heat can be absorbed by a thermostat, 
the ordinary cold crystallization mechanism works [|TJ. So, in the crystallized films of 
semiconductors the undulating surface is developed under the formation of combs, whose 
long axes are perpendicular to the direction of the crystallization front movement || . As 
is known || || , at small number of the crystal embryos this mechanism is realized if the 
nucleation frequency of crystallization centers J is very small. 

Sometimes, explosive crystallization can be initiated by the local heating (for example, 
by laser or electron impulse). Such a scenario takes place in the case of instability appear- 
ance of the interfacial boundary motion due to both the heat exchange with substrate 
and the influence of laser radiation |5], This instability is ensured with a nonlinear 
dependence of the crystal growth velocity u as a function of temperature. Moreover, the 
crystallization front instability can be fluctuational in character that manifests itself in 
the experiments with undulation of crystallized surface || Such a behaviour appears 
in the partially crystallized film or at the incomplete crystallization of amorphous phase. 

Another scenario is observed in the amorphous medium in which as the crystal growth 
velocity u, so the nucleation rate J of embryonic crystals are sufficient large. Indeed, at 
low temperatures, the quantities J and u increase with temperature growth, so that the 
self-heating stimulates the crystallization. Therefore, the increase of a film thickness can 
lead to a situation, when the crystallization heat can not be absorbed by an environment 
that causes the heat instability J7J. As a result, the spontaneous transition to regime 
of the explosive crystallization can be provided by the heating effect. Examples of such 
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amorphous mediums are the amorphous ice layers, the some organic matters |3|, and 
the layers of the germanium amorphous ultradisperse powder with the admixture of the 
crystalline phase || £| . 

Our work is devoted to studying of the explosive crystallization mechanism that is 
provided by such type instability. 

The experimental data in Sec. || show that course of explosive crystallization of ultra- 
disperse amorphous materials is determined by high density of crystalline phase embryos. 
The crystallization phase spreading is similar to percolation cluster formation under con- 
sideration of liquid flow in a random medium || . The formed cluster has the branching 
fractal structure that is characteristic for the thermal conductivity limited aggregation. 
Our approach is based on the assumption that such structure leads to the hierarchical 
picture of the explosive crystallization process. 

In Sec. |3] investigation of realization conditions of this picture as a result of self- 
organization is carried out. This process is fixed by the velocity of crystallization front 
motion, its temperature, and difference of the thermodynamic potentials of amorphous 
and crystalline states. The first of mentioned parameters is connected with the third one 
by means of the positive feedback that is the reason for self-organization. The connection 
between the first and second parameters is due to the negative feedback reflecting the 
Le Chatelier principle. As a result of the interplay between the pointed out factors, the 
stationary state is established at supercritical value of transformation thermal energy, 
where the crystallization front velocity can take anomalous large values. Then, at small 
film thickness, a regime is realized when the crystallization process can be initiated only 
by the external influence type of the laser beam 0. However, with film growth up to 
the critical value the crystallization heat in the film volume becomes to be sufficient for 
spontaneous increase of the front velocity. Such a situation is observed experimentally in 
Ref. 1. 

In Sec. H the direct examination of explosive crystallization is carried out as a self- 
organized criticality process caused by stochastic heat spreading over hierarchical tree 
nodes. The study of effective motion equation shows that, in accordance with the re- 
sult of previous section, the instability is developing in the case, when crystallization 
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heat effect (or externally inputted energy) is above a critical quantity determined by the 
thermometric conductivity. The stationary crystallization heat distribution, defined as a 
solution of corresponding Fokker-Planck equation, allows us to find both the heat cur- 
rent arising during crystallization and the probability of spontaneous crystallization in a 
film with subcritical thickness. According to Sec. [|, the probability increases logarithmi- 
cally as a function of time until the maximum value. On its turn, this value decreases 
monotonically with the growth of thermometric conductivity. 

2 Experimental results 

The experimental study of influence of the crystal inclusions distributed in the amorphous 
phase volume in the kinetics of the explosive crystallization has been carried out with the 
germanium Q. Unfortunately, there are no information in literature about spontaneous 
rise of the explosive crystallization in continuous amorphous thin films of germanium. It 
is known only that the maximum value of formation frequency of crystalline embryos is 
Jmax ~ 10 14 s _1 cm -3 in the supercooled germanium and corresponds to the more high 
temperatures than realized at the explosive crystallization process Therefore, in the 
germanium amorphous films the natural process of crystals nucleation has not enough 
intensity for significant influence in the explosive crystal growth. Let us point out in this 
connection that in amorphous ice layers, where the spontaneous explosive crystallization 
takes place, one has J ma x ~ 10 20 s _1 cm~ 3 . 

In order to intensify the mentioned influence of crystal inclusions, the experiments were 
carried out with the thin films of amorphous nanopowders with admixture of large number 
of smallest crystals having not more 3-10% of total mass. The layers of amorphous powder 
with the characteristic size of particles 3-10 nm were obtained by thermal evaporation 
and following condensation of germanium in the atmosphere of inert gas at the pressures 
10-100 Pa. Changing the evaporation intensity allows us to regulate the part of the 
crystalline particles in amorphous powder. Another peculiarity of our experiments is that 
the substrates absorb significantly smaller heat due to porosity of amorphous films. 

The spontaneous explosive crystallization has been observed in layers of nanopowders 
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with the thickness 0.01-0.1 mm at the substrate temperature 300-400 K. In dependence 
on initial concentration of crystalline phase, the movement velocities of crystallization 
wave have been changed in the range 0.01-0.1 m/s. In contrast to the transverse undula- 
tion being inherent in the usual crystallization mechanism, in our case the front movement 
leads almost always to formation of "twigs" along the movement direction. The charac- 
teristic pictures of explosive crystallization in the powder layers is shown in Fig. 1, where 
the light background corresponds to non-crystallized domains. It is seen that the rise of 
explosive crystallization avalanches occurs from the single centers which act as original 
embryos. The cornerstone of our observation is that the crystalline phase distribution has 
a fractal character being similar to the pattern appeared in diffusion limited aggregation 



3 Determination of explosive crystallization condi- 
tions 

The experimental data show that two mechanism of amorphous material crystallization, 
depending on external conditions and presence of crystalline phase embryos, are possible: 
the slow growth of a cold crystal and the explosive crystallization that is caused by the 
phase transformation heating. According to ||, the transition between these regimes is 
jump-like in character, as a first-order phase transition. We will show below that such a 
transition is caused by the system self-organization due to the positive feedback between 
the heating and growth velocity of crystalline phase. 

To analyze the problem, let us examine the time dependencies of the crystallization 
front velocity u(t), its temperature T(t), and the specific crystallization heat /(£). The 
equations defining these dependencies take into account their dissipative character and 
the positive feedback between quantities u and /, that is the reason for self-organization. 
On the other hand, in order to provide the stability of a system we introduce also the 
negative feedback between u and T. The obtained as a result equations coincide formally 
with the Lorenz system that is the simplest way to describe the self-organization process 



The first of the stated equations has the form 

u = -u/r u + fiT, (3.1) 

where the dot stands for a derivative with respect to time t, fi > is a constant. The first 
term in the right-hand side describes the Debye relaxation during time t u , the second one 
reflects the increase of the crystallization front velocity with the growth of temperature 
difference T = T — at the crystallization front and the thermostat, respectively. In 
the stationary state u = and Eq. (|3.1|) gives linear relationship u = A U T, where the 
constant A u = r u fi is introduced. 

The equation for the rate of quantity T variation has the nonlinear form 

f = -T/r T + g T uf, (3.2) 

where / > is the volume density of difference of the thermodynamic potentials of 
amorphous and crystalline states; tt, gr are positive constants. As in Eq. (|3.1|) , the first 
term in the right-hand side of Eq. (|3.2|) describes the relaxation process of temperatures 
difference T to the stationary value T = 0. It takes place not during the macroscopic time 
r u but during the mesoscopic one tt, so that the important for the future condition tt <C 
t u is satisfied. The second term describes the mentioned positive feedback between the 
crystallization front velocity u and the difference / between the specific thermodynamic 
potentials of the phases that results in the increase of value T and, thus, causes the 
self-organization process. In the stationary case T = Eq. (|372|) takes form 



T = + A T uf, (3.3) 



where the constant At = Qt t t is introduced. According to Eq. (|3.3| ) nonlinear term 
describes the heating of the crystallization front with the growth of crystallization wave 
velocity. 

The kinetic equation for the difference / of specific thermodynamic potentials 

/ = (/o - f)/r f - 9fuT (3.4) 

differs from Eqs. ( |3.1|) , (|3.2| ) as follows: the relaxation of quantity / occurs not to the 
zero but to the finite value /o, representing the energy density inputted in the system 
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(heat effect of transformation); Tf is a corresponding relaxation time. In Eq. ( |3.4j ) the 
negative feedback between the quantities u and T is introduced to imply the decrease 
of thermodynamic transformation factor / with the growth of the crystallization front 
velocity and its temperature (gj > is a corresponding constant). 

Let us study the system of differential equations (|3.1| ), (|3.2| ), (|3.4j ) defining the self- 
consistent behaviour of the quantities u(t),T(t), f(t) which act as the order parameter, 
the conjugate field, and the control parameter, respectively 0. With this aim, we write 
the Lorenz system in the form 

r u ii = -u + A U T, 

r T f = — T + A T uf, (3.5) 
T f f = (/o - /) - A f uT, 

where the constant Af = Tfgj is introduced. As is known |J, the analytic examination is 
possible only provided that hierarchical subordination conditions are satisfied 

T T < T u , Tf < t u , (3.6) 

which mean that in the course of evolution the temperatures difference T and the ther- 
modynamic potential / follow the variation of the crystallization front velocity. As was 
mentioned above, the first of these conditions is always obeyed. Since, on the other hand, 
Tf ~ r T the second inequality (|3~6|) is met also. 



When the values of relaxation times r u , tt, Tf are constant and the conditions ( |3.6| ) 
are obeyed, it is not difficult to see that the system of equations ( |3.5| ) describes the second- 
order transition. However, the cold crystallization mode transforms into the explosive one 
in accordance with the first-order mechanism. To avoid this discrepancy, let us use the 



simplest approximation |T0| 



- = - (1 + W|, (3.7) 

characterized by the positive constants tq, k, and u T . Moreover, it is convenient to intro- 
duce the scales of quantities u, T, and /: 

u m = {A T Af)- l '\ T m ee u m /A u = A- u \A T Af)- l '\ f c = {A U A T )~\ (3.8) 
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where A u = tq/x. Then, the system ( |3.5| ) assumes the simplest form 



t u = —u 



I + n{l + u 2 / a 2 Y l \ + T, (3.9) 
TT f = — T + uf, (3.10) 
T f f=(fo-f)-uT, (3.11) 

where the parameter a = u T /u m is introduced. 

Taking into consideration conditions (|3.6|) , the left-hand sides of Eqs. (|3.10|) , (|3.11| ) 



can be set equal zero because of the small relaxation times tt, tj. As a result we obtain 
the expressions for the temperatures difference T and the thermodynamic potential / in 
terms of the velocity u of crystallization front: 

T = f u(l + u 2 )-\ (3.12) 
/ = f (l + u 2 y\ (3.13) 



At u <C 1, the dependence ( |3.12|) has the linear form characterized by susceptibility 
du/dT = f^ 1 . At u — 1 function T{u) becomes saturated, and at it > 1 decreasing 
dependence is realized that has no physical meaning. This implies that the constant u m 
defined in ( |3.8| ) has the physical meaning of the maximum value of the crystallization front 
velocity. According to Eq. ( |3.13D the difference / of specific thermodynamic potentials of 
different phases decreases monotonically with the growth of velocity u from the value fo 
at it = to the / /2 at u = 1. Obviously, this decrease is caused by the negative feedback 



in Eq. ( |3.11| ) that is the reflection of Le Chatelier principle for examined problem. On the 



other hand, the positive feedback between the velocity u and the thermodynamic factor 
/ in Eq. (|3.10| ) is the reason for transition to the explosive crystallization mechanism that 
leads to the growth of T due to crystallization front heating. However, in accordance 
with Eq. ( p. 11| ), the latter results in decrease of / to be a consequence of self-organization 



process. 

Within the framework of the adiabatic approximation tt, 7/ to> the Lorenz system 
( |3.9|) - (|3.11|) is reduced to the Landau-Khalatnikov equation 



r u = -dV/du. (3.14) 
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Its form is determined by the effective potential 



,/= 2 



u 2 -9\n (\W) 



na , 
H In 



1+ (u/af 



(3.15) 



where 6 = f /Ah, Ah = (toTt^^)" 1 is the scale defining specific crystallization heat, 
the quantity V is measured in units of u 2 m . For small values of 9 the curve of the V vs u 
dependence has a monotonically increasing shape with its minimum at point u = that 
corresponds to the cold crystallization mechanism. At 9 = 6®, where 



9° c = l + a 2 (K - 1) + 2a Jk (1 - a 2 ) 



(3.16) 



a plateau appears, which for 9 > 9® is transformed into a minimum corresponding to 
the velocity u e ^ and a maximum at point u m separating a minima which meet to the 
values u = and u = u e . When the parameter 9 increases still further, the minimum at 
point u = u e is lowered and the height of barrier at u = u m decreases, vanishing at the 
critical value 



9 C = 1 + K. 



(3.17) 



The steady-state values of the crystallization front velocity have the form 
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As is shown in Fig. 2a, if the system's energy increases slowly, the jump from zero to a/2moo 
is observed at point 9 = 9 C and then the value u e increases smoothly. If the parameter 
9 goes downward quasistatically, the crystallization front velocity u e smoothly decreases 
up to the point, where 9 = 9° c and u e = u o, and then jump-like goes down to zero. The 
hysteresis of such type takes place only at the presence of energy barrier inherent in a 
first-order transition and appears if only the parameter a = u T /u m is smaller than unity. 

The key point of studied transition is that the stationary value of the thermodynamic 
transition factor 

(1 + u 2 ) - y/il + v&f-eil-a?) 



fe 



1 - a 2 



(3.19) 
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equals the thermal energy 9 in the < 9 < 9° c interval (Fig. 2b). At 9 > 9® this quantity 
decreases smoothly from the value 

f m = l + a[K/{l-a 2 )] 1 ' 2 (3.20) 

at 9 = 9° c to 1 at 9 ->• oo. 

Under quasistatic growth of parameter 9 from to 9 C the stationary value of transfor- 
mation factor increases linearly being in the same interval. After jump-down at 9 = 9 C the 
quantity f e decreases smoothly according to dependence ( |3.19| ). Under reverse decrease 



of 9 the quantity f e undergoes the jump at point from the value f m up to the 9®. Since, 
in the important range of values of the parameters a and k limited by n m in — & 2 / (1 — « 2 ) ; 
the maximum value f m of specific transition energy is smaller than the minimum value 
9® of the heat density, the stationary value f e of the specific difference of thermodynamic 
potentials of amorphous and crystalline states is always smaller than heat density 9. 

The above analysis shows that the effective potential V(u) has the barrier separating 
the cold and the explosive crystallization modes. As heat density 9 becomes greater than 
the critical value 9 C this barrier disappears. Thus, at 9 < 9 C the transition to the explosive 
crystallization mechanism requires the penetrating of energy barrier and at opposite case 
it realizes spontaneously. The first of appointed situations takes place in the case when 
the explosive crystallization is initiated by an external beam (see Fig. la). With the 
increase of the coating thickness the crystallization heat can not be absorbed by substrate 
and parameter 9 increases. This leads to that the value 9 C (at which the function V(u) 
loses barrier) is reached at the critical film thickness and the system transforms into the 
explosive crystallization regime spontaneously (Fig. lb). 

4 Description of explosive crystallization as a self- 
organized criticality process 

In recent years considerable study has been given to the conception of self-organized criti- 
cality (SOC) representing the natural development of critical phenomena picture [IT, O. 



The basic distinction of the SOC from the phase transition is that SOC process realizes 
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spontaneously, whereas the phase transition goes on under the external influence (for 
instance, the temperature variation). In this section we will expound the quantitative 
picture within the framework of which the explosive crystallization will be represented as 
SOC. The basis for the such representation is that the SOC process consists of the hier- 
archically subordinated sequence of elementary phase transitions which are usually called 
avalanches jl3" l. Their hierarchical subordination manifests itself in that the avalanches of 



upper level are formed after their formation has been finished on the lower one. Then this 
process recurs on the more upper levels - right up to the global avalanche formation on 
the top of hierarchical tree. The hierarchical nature of explosive crystallization process is 
discovered obviously in the microscopic photographs of crystallization pattern in Fig. 1, 
where it appeares as the tree-like fractal structure. In accordance with |L4| the hierarchi- 
cal tree represents the geometrical shape of ultrametric space in that system's states are 
realized. Thus, for description of explosive crystallization the geometrical picture of the 
nodes distribution over hierarchical tree levels is necessary to represent, at first . 

Let maximum number of nodes N be on the bottom hierarchical level corresponding 
to the distance in the ultrametric space s = 0. This level meets the elementary avalanches 
whose number coincides with N. There is the only node on the top level (s = So > 1) 
corresponding to the global avalanche. The problem is to find the dependence N(s) that 
define the distribution of tree node number over hierarchical levels. 

At first, we examine the basic types of the trees (Fig. 3): regular tree with integer 
branching ratio j, regular Fibbonachi tree with fractional j = r « 1.618, degenerate tree 
with the only branching node per level and the tree of our primary concern - irregular 
tree. Let k be the numbering index for the levels, so that k increases from the top level 
to the bottom one. The variable 

s = s Q -k (4.1) 

then defines the distance in the ultrametric space |14|, |i~6| . Geometrically, objects of this 
space correspond to the nodes of the bottom level {k = s ) of a Cayley tree. Since the 
distance between the nodes is defined by the number of steps to a common ancestor, the 
distance is eventually the level number (|4.1|), counted from the bottom. 
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As it can be seen in Fig. 3a, in the simplest case of regular tree with integer branching 
ratio j the number of avalanches N k = j k exponentially decays to zero with the distance 
s between them: 



N(s) = JVexp(-slnj), N = j s °. (4.2) 

In Eq. Q4.2|) the equality (|4.1|) is used and the avalanche number N is related to the 
total number of levels sq. For the Fibbonachi tree (see Fig. 3b), where N k = vr k , v ~ 
1.171, r « 1.618, we have 

N(s) = iVexp(-slnr), N = ut s °. (4.3) 



When Eq. (|4.3|) is compared with Eq. ( |4.2| ), it is clear that the exponential decay remains 
unaltered in the case of fractional branching ratio and characterizes the regularity of tree. 

For the degenerate tree (see Fig. 3c) iV fe = (J — l)k + 1 and Eq. (|4.1|) provides the 
following linear dependence 

N(s) =N-(j- l) a , N=(j- l)s + 1. (4.4) 

It can be shown that in the case of irregular tree, displayed in Fig. 3d, the power law 
dependence is realized: 

N k = k a , a> 1. (4.5) 



The latter can be regarded as an intermediate case between the exponential Eqs. (|4.2|) , 
( f4.3| ) and linear Eq. ( |4.4| ) obtained for the limiting cases of regular and degenerate trees, 
respectively. Formally, the approximation (|4.5|) means that a function N(x) defined on 
the self-similar set of hierarchically subordinated avalanches is homogeneous, N(kx) = 
k a N(x). It is convenient to rewrite Eq. ( |4.5| ) in terms of the distance: 

N k = N(l - s/s ) a , N = s a , a>l. (4.6) 

At given value of crystallization thermal effect Q k the heat current density between 
different levels k is expressed by the generalized Onsager equality 

3k = -x(Qk)^r- (4-7) 
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Here, within the multiplicative noise representation [TJJ the effective thermometric con- 
ductivity coefficient 

x(Q) = xQ p (4.8) 

is defined by the constant \ > and the exponent (3. The cornerstone of our approach is 
that total heat current at given level is independent of the hierarchical level: 

j k N k = const = J. (4.9) 



Substitution of Eqs. (|4.6p - (|4.8| ) into Eq. (|4.9| ) gives the expression for the crystallization 
heat effect: 

Q k = Qk~\ 6=(a-l)/(l + /3) >0 (4.10) 
normalized by the maximum value Q = Q k =i- Inserting Eq. ( |4.1|) , we get the dependence 

Q(s)=q(l-s/s )- b , (4.11) 
where the heat effect at the bottom level s = is 

q = Qs b = QN- b ' a . (4.12) 



In general case, the condition of current conservation ( |4.9[ ) is not satisfied and with 
accounting Eqs. ( |4.10| ), ( [1.12| ) we assume the scaling relation 



Qk = N h ' a k- b q kl (4.13) 

where q k is a slowly varying function to be determined. According to Eqs. ( |4.7] )-( ]4~9| ) it 
obeys the Landau-Khalatnikov equation 



dx _dV 

d.K dx 



(4.14) 
where one denotes 

« = In k b , x = q k /q c , q l c +P = (J/b X ) N'^ a , (4.15) 
and the effective potential is 
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As indicated in Fig. 4, this potential reaches its maximum value Vq = (1 +/3)/2(l — (3) 
at x = 1 and decreases indefinitely at x > 1. So, in order to initiate the process of 
explosive crystallization, a low intensity avalanche with q < q c at the bottom level needs 
to penetrate the barrier Vq. For study of this process we proceed with stochastic Langevin 
equation with a white noise (cf. Eq. (|4 . 1 4j) ) 

dx dV ^ , A 

(0 = 0, (C(«)CM> = 2 X 5(«- K '), (4.18) 



where the noise intensity x is reduced to the thermometric conductivity in Eq. 

The solutions of this equation are distributed in the ultrametric space according to 
the function w(n,x) that obeys the Fokker-Planck equation 



dw dj dV dw 

<7K OX OX OX 

Since there is no current at the equilibrium state (j = 0), the distribution function 

w (x) oc exp (-V(x)/x) (4.20) 



is dictated by the potential (|4.16|) . In the case of non-equilibrium steady state the prob- 



ability density w does not depend on the hierarchical level variable k and the current j 
being constant, in compliance with conservation law ( |4.9| ), can take a nonzero value. In 
accordance with Eq. (|4 . 1 9|) the stationary w(x) and the equilibrium Wq(x) distributions 



are connected by the equation |L9 
w(q) j f dx 



w (q) X / w (x) 



(4.21) 



where the boundary condition w — > at q —>■ oo is taken into account. 

Given the heat effect q Eq. ( [4.21 ) allows the current j to be found. In trying to do it, 



special consideration should be given to the fact that the heat q is bounded from below, 
q > G [JT2"| . The appearance of the gap G is the feature inherent in hierarchical ensemble 



of crystallization centers. Indeed, after merging of them within a hierarchical cluster of 
the size s g , all s, such that s < s g , are appeared to be dropped out the consideration as 
well as low heat effects with q < q(s g ) = G (see Fig. 3). The expression for the current j 
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then can be derived from Eq. ( |4.21| ) with the second boundary condition w{G) = wo(G). 
The result reads 



j = 2 X W 
where the factor 



-i 

(4.22) 



W oc exp(-V /x), V ° = \ \z4 ( 4 - 23 ) 

gives the probability that fluctuation will surmount the barrier Vq of the potential ( |4.16|) . 
Eq. (|4.22| ) shows that in the case of small gap, G <C q c , the current j has the value W\ 
and it is doubled under G = q c . It can be understood if we picture the effect of the 
gap as a mirror that reflects diffusing particles at the point q = G: if G <C q c a particle 
penetrating the barrier can move along both directions, but in the case of G = q c the 
mirror is placed at the point corresponding to the top of the barrier and all particles go 
down the side where the quantity q grows indefinitely. 

Given the current j the stationary distribution function w(x) is defined by Eq. ( |4.21| ), 
according to which, w(x) ~ wo(x) in the subcritical region q < q c , while in the supercritical 
range q ^> q c we have w (x) ^> w(x) due to indefinite increase of w (x). As far as the 
stationary distribution is concerned, it can be derived from the current definition (|4.19|) , 
where the last diffusion term is negligible for supercritical heats: j ~ — (dV/dx)w. The 
result is that the probability w(q) remains almost unaltered, w(q) ~ w(q c ), in the range 
from q c up to a boundary value q g and w(q) ~ at q > q g PUL The growth of q g is 
governed by the equation 

^ = x^- (4-24) 
dn q 2 g 

Since the above picture is essentially statistical, it enables the critical heat effect q c for 
the transition point to be found. Indeed, when the definition of the macroscopic current 
J in Eq. ( |4.9| ) is compared to that of the microscopic current j in Eq. ( [4.19[ ) , it is apparent 
that they differ from one another only by the factor N^ 1 ^ 0, = dependent on the 
total number of embryos N (see Eq. (|4.6| )). On this basis, the last expression of Eq. ( |4.15| ) 
and Eq. ( |4.22| ) at G = 0, x ^ 1 gi ye the desired result: 



g c = Qexpf- (1 f^ 1 ), (4.25) 



2X 
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where the pre-exponent factor Q determines the probability of the barrier penetrating 
and cannot be calculated within the framework of the presented approach. Eq. (|4.25| ) 
predicts the slow growth of the critical thermal effect q c of embryo with the thermometric 
conductivity coefficient x- 



5 Time dependence of crystal formation probability 

Since the ensemble of hierarchically subordinated crystallization centers represents a self- 
similar set, the probability distribution of embryos P(Q, s) in the course of SOC process 



is a homogeneous function of s p~2 



P(Q,s) = s- T w(q), (5.1) 

where w(q) is the stationary distribution of embryos and r is the positive exponent. 
Physically, Eq. (|5.1|) implies that the heat effect Q, being measured by the scale N b / a , 
equals the heat effect of an embryo formation q in accordance with Eq. Q4.12 ). 

In this section we are aimed to define the probability of hierarchical crystallization 
leading to the formation of a fractal cluster (see Fig. 1). As it has been clarified in 
Sec. |], this process can be conceived of as diffusion in ultrametric space that makes 



the distribution (|5.1|) mounted. In order to find the conditional probability V(t) that 
no crystallization will begin at time t one has to integrate over s the distribution (|5.1| ) 
weighted with the function 

p s (t) = exp (-t/t(s)), t(s) = t exp [Q{s)/ X \ (5.2) 

descriptive of Debye relaxation with the time t(s) governed by the barrier height Q(s) (see 
Eq. ( |4.11| )) and a microscopic time to determined below. By using the steepest descent 
method, it is not difficult to derive the late time (t 3> t e /) asymptotic formula 



V(t) 



q 



r/b 



1 



-I/ft' 



t. 



ef 



(5.3) 



.1 tef, 

Eq. ( |5.3j ) has been obtained by assuming that the condition 1 s m < so is met, where 
s m denotes the location of the maximum of integrand and obeys the equation 



XT (I 

bq 



x 



x 



,1+6 



t 

T 



exp 



X 



l-x) 



X = 



so 



(5.4) 
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Taking into consideration the scaling relation for the number of hierarchical levels sq, 



which is the cut-off parameter [12 



s oc (q c - q) 1/(7 , o > (5.5) 
we readily come to the conclusion that the condition is satisfied provided 

g-g c <g, t > t e/ exp ([q c lxY l,b - i) • (5.6) 

Clearly, from Eq. ( |5.6j ) the heat effect q in Eqs. (|5.3|), ( |5.4| ) can be replaced by q c . Note that 
in accord with Eq. ( |5.3| ) the crystallization probability V(t) = 1 — V(t) logarithmically 
increases in time up to the value V = 1 — {q c /xY^ h - The condition V > is satisfied 
if in Eq. ( f4.25|) factor Q = (e/2)(l — /3) _1 and thermometric conductivity coefficient is 
restricted by the maximum value Xo — (1/2) (1 — j3)~ l . 



6 Discussion 

In accordance with above approach, given in Sec. [^, the transition from the cold crystal- 
lization mode of amorphous material to the explosive crystallization mechanism represents 
the self-organization process realized as a first-order transition. The crystallization front 
velocity u represents the order parameter, the temperatures difference T at the crystal- 
lization front and thermostat acts as the conjugate field, and the difference / of specific 
thermodynamic potentials of amorphous and crystalline states is the control parameter. 
Equations (|3.9|)-( |3~TTT| ) are derived assuming the degrees of freedom u, T and / to be dis- 
sipative. In addition, the positive feedback between u and / is taken into consideration 
as the reason behind the self-organization, whereas the negative feedback between u and 
T is a manifestation of Le Chatelier principle. The system is driven by the parameter 
6, whose value represents the thermal energy of crystallization (or externally inputted 
energy). When 6 is above the critical value ( |3.16| ), the effective potential (|3.15|) assumes 
additional minimum value — q < at point u = u e and maximum one U at u = u m (sta- 
tionary values of crystallization front velocity u and thermodynamic factor / are given 
by Eqs. ( |3.18|) , ( |3.19| )). The explosive crystallization process is preferable in potential 
provided the minimum effective potential becomes negative (q > 0). The height of the 
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energy barrier U of the effective potential ( |3.15| ) defines characteristic time in the last 
Eq. ( |5.2|) for hot crystallization center to be formed 



t « T D exp(U/A), 



(6-1) 



where A is the variance of 9 and ~ 10~ 12 s is the Debye time. 

The nucleating crystals form statistical ensemble of hierarchically subordinated ob- 
jects, characterized by heat q and distances s in ultrametric space (crystalline cluster 
size |12[]). Since the crystallization represents the effective diffusion over hierarchical tree 
nodes, then, similar to Brownian particle with coordinate q at time s, the ensemble can 
be described by Langevin equation ( |4.17|) subjected to the noise Eq. ( f4.18|) with \ being 
the effective diffusion coefficient (thermometric conductivity) and corresponding Fokker- 
Planck equation (}4.19|) . The stationary heat distribution and the steady-state current are 
given by Eqs. ( |4.21| ), ( |4.22j ). The condition of current conservation Eq. ( f4.9|) yields the 
heat distribution ( f4.11| ) in the ultrametric space. The ensemble of embryos, being weakly 
dependent on s, is governed by the effective potential (|4.16|) that reaches its maximum at 
the critical heat ( |4.25| ) (see Fig. 4). So, the explosive crystallization requires supercritical 
heat effect, q > q c , to surmount the barrier Vq with the characteristic time (cf. Eq. ( f4.23| )) 



T^t exp (V /x). 



(6.2) 



This picture bears some resemblance with the formation process of supercritical embryo in 
the theory of a first-order phase transitions |B1. In the course of phase transformation the 
next stage is the growth of the embryo due to the diffusion increase of the heat effect Q(s) 
in ultrametric space. As a result, we have the logarithmically slow large time asymptotic 
for the probability of the total cluster formation: 

-1/61 - T 



V(t) 



V 



i-Wh,'- r 



t. 



where time t is counted from the instant T, Eq. 
that no crystallization will occur: 

r/b 



exp 



t ef = (r/b)V 1/T t , (6.3) 
and V is the minimum probability 

(6.4) 
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From Eq. ( |6.4| ) the probability is determined by the ratio of the noise intensity x ( see 
Eq. ( (4.18|) ) and its maximum value xo = (1/2) (1 — /3)~~ 1 - The key point is that the 
maximum probability V = 1 — V of crystallization is completely suppressed under great 
thermometric conductivity coefficient % (Fig- 5). 
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FIGURE CAPTIONS 



Fig. 1. The patterns which arise in the layer of amorphous powder of Ge at the ex- 
plosive crystallization: (a) from single center; (b) from several centers which arise 
spontaneously. 

Fig. 2. (a) The dependence of the steady-state velocities of crystallization front on the 
transition heat (the solid curve corresponds to the stable state u e , the dotted curve 
meets the unstable one u m ). (b) The dependence of the stationary difference f e of 
specific thermodynamic phase potentials on the transition heat. The arrows indicate 
the hysteresis loop. 

Fig. 3. Different types of hierarchical trees (the level number is indicated at left, cor- 
responding number of nodes - at right): a) regular tree with j = 2; b) Fibbonachi 
tree; c) degenerate tree with j = 3; d) irregular tree. 



Fig. 4. The form of the effective potential ( ^.16| ) at j3 = 0.2. 



Fig. 5. The dependence of the maximum probability V of crystallization on the thermo- 
metric conductivity \- 
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